Size-dependent thermal conductivity of nanoscale semiconducting systems 
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We study the size dependence of thermal conductivity in nanoscale semiconducting systems. An analytical 
formula including the surface scattering and the size confinement effects of phonon transport is derived. The 
theoretical formula gives good agreements with the existing experimental data for Si and GaAs nanowires and 
thin films. 
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It is of great importance from both fundamental and ap- 
plication point of views to understand heat conduction in 
nanoscale materials. On one hand, many fundamental ques- 
tions remain open, such as whether the thermodynamical laws 
still apply to the nanoscale materials which are characterized 
by finite number of atoms/molecules, while the thermodynam- 
ics is established on the fact that the system contains infinite 
large number of molecules. On the other hand, a continuous 
miniaturization and increase of running speed of semiconduc- 
tor devices will give rise to redundant heat that will deteriorate 
the devices' efficiency.' How to dissipate the heat from these 
devices becomes a crucial issue in heat management. 

Generally, two approaches are used in studying heat con- 
duction in such low-dimensional systems of finite size. At 
first, in order to understand the physical mechanism of heat 
conduction in such systems, many one-dimensional (ID) lat- 
tice models with and without on-site (pinning) potential are 
used such as the Fermi-Pasta-Ulam (FPU) mode!^ without on- 
site potential, and the Frenkel Kontoroval (FK) model and the 
(f)'^ model with on-site potential.^ More recently, a polymer 
model with transverse motion has been introduced and studied 
both analytically and numericallyii All these models nicely 
demonstrate the role of anharmonicity and the role of on-site 
potential in ID heat conduction. They are useful in helping 
us understand the fundamental law of heat conduction, the 
Fourier law. In the system without on-site potential, a size- 
dependent thermal conductivity has been observed,^ which is 
attributed to a superdiffusive motion of phonons.^ However, 
such models are too simple to be used for modelling heat con- 
duction in realistic nanostructures, because the heat conduc- 
tion in such models is restricted to ID or quasi- ID, the scat- 
tering of phonons from surfaces is completely neglected. 

Another popular approach employs the Boltzmann trans- 
port equation, which can explain the bulk transport phenom- 
ena well. However, as we know that for nanostructures, 
both the effects of phonon confinement and non-equilibrium 
phonon distribution due to boundary scattering are impor- 
tant. Indeed, recent theoretical and experimental investigation 
demonstrates a size-dependent behavior of lattice vibration of 



nanostructures.^ Experimental measurements also show that 
thermal conductivity of semiconducting nanowires and thin 
films is smaller than the corresponding bulk value. How- 
ever, a general satisfactory theory to give a prediction agreed 
with existing experimental results at room temperature is not 
available yet. 

In this paper, we propose a phenomenological theory for 
the size dependence of thermal conductivity by taking into 
account the intrinsic size effect of phonon velocity, mean free 
path and the surface scattering effect. An explicit analytic ex- 
pression is obtained, in which all parameters have clear phys- 
ical meaning, and good agreements with the existing experi- 
ments are found. The different roles of two basic effects in the 
different size range are illustrated. 

We start with the well-known kinetic formula of thermal 
conductivity k for bulk dielectric materialspi 



K = —Cvl, 
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where c is the specific heat, v the average phonon velocity and 
I the mean free path (MFP). We assume that phonons predom- 
inate the thermal conduction, which is true for semiconduc- 
tors or insulators discussed here. For thermal conductivity of 
nanomaterials, we first consider the size dependence of v and 
I. The specific heat is assumed to be constant and the room 
temperature is applied in this paper. The average phonon ve- 
locity is proportional to the characteristic Debye temperature 
of crystals|i2ii^ 
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with the Planck constant h, the Boltzmann constant ks, the 
Avogadro constant Na, and the molar volume V . The system 
considered is assumed to be isotropic. 

In the following discussion, we suppose Eq. (|3 is valid 
for the corresponding nanoscale crystals. Let L be the size 
of nanostructures, such as the diameter of nanowires or the 
thickness of thin films, the L ^ oo limit is denoted by using 
the subscript h which means the corresponding bulk limit, the 
size dependence of the phonon velocity is equal to that of the 
Debye temperature, we thus have, vl/v\j = Ol/O^. 

From the Lindemann's proposition, we may get the rela- 
tionship between the melting point and the Debye temperature 
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of crystals. In 1910, Lindemann proposed a melting criterion, 
known to be valid for small particles, stating that a crystal 
will melt when the root mean square displacement (MSD) <j 
of atoms in the crystal exceeds a certain fraction of the inter- 
atomic distance.'^ Combining with the Einstein specific heat 
theory, the square of the characteristic temperature is propor- 
tional to the melting point T™ of crystals, and the modern 
form of this relation for the Debye temperature isiiiii 

/ rpm \ 1/2 

® = ™- \-Mvm ) 

with the molecular mass M. According to the same relation 
for nanocrystals, 

Ql/Ql = T^/T^. (4) 

Combining the Lindemann melting formula and the Debye 
expression for thermal conductivity, eliminating the Debye 
temperature and the sound velocity, it can be deduced that the 
MFP is proportional to the melting point at a given absolute 
temperature T, I cx 20r'"a/(7^T) with the lattice constant a 
and the Gruneisen constant 7." '^ Similarly, the same relation 
is assumed for nanocrystals, we have 

h/lb = T^/T^. (5) 

From Eqs. Q, @ and (|5Jl, we obtain an expression for 
the size-dependent thermal conductivity of nanosemiconduc- 
tors [y™/y™]'^/2 Moreover, based on the fact of the size- 
dependent atomic thermal vibrations, the relations among the 
atomic MSD cr^, the Debye temperature and the melting point, 
the size-dependent melting temperature function of nanocrys- 
tals has been modelled and validated as followsi^ 

Tr/rr = {cr,/aLf = exp ( ^j^o^-j ) , (6) 

where a{= a^j erf) is a material constant with erf and erf cor- 
responding to the MSD of surface atoms of a crystal and that 
of atoms within the crystal, respectively. For the free-standing 
nanocrystals, a = 2S'„/(3i?) + 1 > 1 with the bulk vibra- 
tional entropy Sy of melting and the ideal gas constant i?"^ 
based on the Mott's expression for the vibrational entropy and 
its relation with the melting point. '^ '^ Lq is a critical size at 
which almost all atoms of a crystal are located on its sur- 
face, Lq — 2(3 — d)w with the atomic/molecular diameter 
w and the dimension d = 0, 1, 2 for nanoparticles, nanowires 
and thin films, respective ly3 Note that the two basic assump- 
tions in the model of the size-dependent atomic thermal vi- 
brations of nanocrystals are: (1) although er^^ and er^^ are 
considered to be size-dependent (the phonon softening is con- 
sidered to occur not only on the surface but also in the in- 
terior for small size crystals), a, the ratio between them, is 
taken approximately as a size-independent value; (2) the co- 
operative coupling between the surface region and the interior 
region is considered phenomenologically by taking the vari- 
ation of er^ to be dependent on the value of itself, er^ is the 
average MSD over the crystal with the respective weight of 



TABLE I: Parameters used in Eq. Q for predictions in this paper. 
a = 2Sy/{3R) + l,Sv^ Sm - -R (Ref. 22) and = -ffm/T™ 
with the bulk melting entropy Sm, enthalpy Hm and temperature 
T™. I/O = 4w, 2-10 for nanowires and thin films, respectively. 





i7„(KJmori) 


r,„(K) 


w (nm) 


lo (nm) 


Si 


50.55^^ 


1685^^ 


0.3368-"* 


41-1 


GaAs 


12C'* 




0.245" 


5.^ 



(Js^ and erj^J^ With the above consideration, a change in cr^ 
can be given by cr^(x + dx) — cr^(x) = (a — l)a^{x)dx, 
where the surface-volume ratio x — Lq/{L — Lq). By inte- 
grating the above equation, the size-dependent MSD was ob- 
tained as shown in Eq. (|6}. According to er^ cx T/Q^ in the 
high temperature approximation {T > 6/2)i^i2£ and Eq. @, 
c^bM = ©i/Qfc = T^/Tl'\ Equation © indicates that 
the melting point and thus the thermal conductivity decreases 
with reducing size L of nanocrystals, it is valid for crystals of 
L > 2Lo with 2Lo representing a characteristic length scale 
for the crystallinity.'» 

The above discussion on the thermal conductivity of nanos- 
tructures, based on the effective bulk formula, includes the 
intrinsic size effect of the phonon velocity and the MFP. The 
phonon-phonon interaction increases with size reduction due 
to the confinement, which causes the increase of thermal re- 
sistance and the decrease of heat conduction. On the other 
hand, as the size decreases and the surface-volume ratio in- 
creases, the large surface/interface scattering, corresponding 
to certain boundary conditions in the linearized Boltzmann 
equation, has great influence on the transport. Considering 
the nonequilibrium phonon distribution due to boundary scat- 
tering, the effect of the surface roughness with the boundary 
scattering shows an exponential suppression in the distribu- 
tion and the conduction."-^' Therefore, a term p exp{—lo / L) 
is added to correct the bulk formula, where p is a factor reflect- 
ing the surface roughness, Iq is the phonon MFP in the Debye 
model at room temperature and assumed to be a constant here 
since we have considered the size effect in the above discus- 
sion. Iq /L corresponds to the Knudsen number of the phonon 
Knudsen flow induced by the interface scattering.^' Finally, 
the size-dependent thermal conductivity is obtained as 

This is the central result of this paper. It gives a quanti- 
tative prediction for the size-dependent thermal conductivity. 
< p < 1. The larger value of p corresponds to the smaller 
roughness, i.e. the smoother surface, thus the more proba- 
bility of specular scattering, vice versa, the smaller p corre- 
sponds to the more probability of diffusive scattering. The 
smaller L corresponds to the relatively larger surface rough- 
ness usually and thus the smaller p. p depends greatly on fab- 
rication precision of nanostructures, its physical meaning will 
be discussed further later on. When system size deceases, the 
value of exp{—lf)/L) decreases, which reflects the increase 
of the interface scattering and thus the weakened conduction. 
From the point of view of the bulk approach, there are two 
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FIG. 1: Size-dependent dimensionless thermal conductivity for Si 
films in terms of Eq. 0, «:l/k6 versus system size L for p = 0.1, 
0.7 and 1, respectively. The other parameters in Eq. Q are given in 
Table I. 
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FIG. 2; Size-dependent thermal conductivity of Si nanowires and 
thin films at 300 K. The symbols are experimental results, the up- 
triangles, circles and down-triangles for thin films (in-plane) cited 
from Refs.9, 10 and 29, respectively, the squares for nanowires (along 
the axis),* The curves are predictions from Eq. 0, the upper one for 
thin films with p — 0.78, the lower one for nanowires with p — 0.4. 



asymptotic limits to be satisfied by Eq. Q: L — > oo and 
p 1, kl L ^ Lq or p — *■ 0, — > 0. In Fig. 

we show kl I Kb versus system size L for different roughness 
parameter p. The thermal conductivity spans a wide range, in 
which our prediction at p = 0.7 for L < 100 nm is consistent 
with the previous theoretical prediction.^' The intrinsic size 
effect decreases the conductivity obviously when the size is 
smaller than about 100 nm. The surface effect related to p can 
explain why the conductivity is low even in the large size in 
some case. 

Figure |2] shows the size-dependent thermal conductivity of 
Si nanowires and thin films from the experiments and the com- 
parison with our theoretical predictions from Eq. (0 with suit- 



able parameters p. The predictions are in good agreements 
with the experimental results. For carefully prepared Si films 
with thickness of 70-250 nm,** p = 0.78 corresponds with 
the experimental results, this large value of p implies small 
diffusive scattering contribution due to the smooth surface re- 
lated with the careful preparation of the nanostructures. For Si 
nanowires with diameter of 20-100 nm,^ p = 0.4, reflecting 
the larger surface roughness of the nanowires in the sample 
fabrication. 

Now we turn to the underlying physical meaning of the 
roughness parameter p. For the Si films with the thickness 
L = 100 nm, if the surface roughness rj, the mean root square 
deviation of height of the surface from the reference plane, is 
2.2 nm, which may be the case according to the precise fabri- 
cation technology of monocrystalline Si layers with thickness 
dispersion smaller than 4 nm^ffi p — 0.78 can be obtained 
by a relation p — 1 — lOrj/L; according to this relation, if 
p — 0.4, f] is 1.32 nm for L = 22 nm, which agrees with 
the experimental observation for the Si nanowires well,^ thus, 
this fitted formula is valid and can be used to determine p by 
measuring or estimating 77 of samples. The larger 77 for a given 
L corresponds to the smaller p, i.e., the rougher surface corre- 
sponds to the higher probability of diffusive scattering. There- 
fore, carefully fabricated structures with smoother surface will 
bring forth a larger conductivity under the same other condi- 
tions. 

Figure |2] also shows that p = 0.78 is applicable for the 
films spanning the size range of 20-100 nm, 70-250 nm and 
400-1600 nm, since the decreased rj can be achieved with de- 
creasing L. Note that the surface roughness may be different 
for samples of different size even in the same experiment, but 
it can be controlled to be fewer than tens of nanometers gen- 
erally in modern fabrication.^" 

Figure|3]shows the thermal conductivity of GaAs nanowires 
and thin films. The theoretical predictions with corresponding 
p also show agreements with the experimental results from 
different groups. Similarly to Si nanowires, p = 0.07 is also 
small for GaAs nanowires,^ kl < 0.1/t;,. It seems that for 
nanowires, the thermal conduction is much worse compared 
to that of the corresponding films, which may be attributed 
to the restrained dimension, the surface fabrication precision 
and thus the increased surface scattering. For the superlattices 
with thickness of 140 nm,-" p = 0.98, kl ~ 0.9Kb is very 
high. For the superlattices of 10-100 nm,-'^ p = 0.75 and kl 
is lower. The smaller L corresponds to the relatively larger 
surface roughness and the smaller p, thus the lower conduc- 
tivity. For the superlattices of 5-45 nm^p — 0.35, kl corre- 
sponds to the cross-plane conductivity measured and is even 
lower. Note that this low conductivity results also from the 
scattering contributions of the multi-interface corresponding 
to the increased interface roughness 77, which has the weaker 
effect on the in-plane conductivity. When the size decreases to 
about 5 nm, the structure and energy state of these nanosolids 
may be different from that of the corresponding crystals, the 
small size and large surface effect may not be sufficient to de- 
scribe the physical properties of such small size structures and 
the quantum effect may be dominant, the present model is not 
suitable for that case. Some molecular dynamic calculation 
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FIG. 3: Size-dependent thermal conductivity of GaAs nanowires and 
thin films. The star is experimental result of nanowire with diame- 
ter about 180 nm (corresponding to the rectangle cross-section with 
the same area) at 40 KJ considering the similar conductivity dif- 
ference between the nanowire and the bulk extrapolated at 300 K, 
the curve is prediction from Eq. Q with p — 0.07. The dia- 
mond is experimental result of GaAs/AlAs superlattice with thick- 
ness 140 nm,"" assumed as the conductivity of a single GaAs film 
with the same thickness considering the similar structures between 
GaAs and AlAs (the following is the same), the curve for prediction 
with p — 0.98. The squares are experimental results of GaAs/ AlAs 
superlattices with thickness 10-100 nm (in-plane),''' the curve for 
prediction with p — 0.75. The circles are experimental results of 
GaAs/ AlAs superlattices with thickness 5-45 nm (cross-plane),'^ the 
curve for prediction with p = 0.35. 



shows the minimum with the size in the cross-plane thermal 
conductivity for the thinner superlattices without roughness, 
which is also theoretically attributed to a crossover between 
the particle transport and the wave transporti^ 



In summary, we have derived a quantitative formula for 
the size-dependent thermal conductivity of nanoscale semi- 
conducting systems by taking into account the intrinsic size 
effect and by correcting the bulk expression. The formula 
relates the thermal conductivity with the system size by the 
surface roughness parameter, the Knudsen number, the crys- 
tallinity length scale and the atomic vibration parameter. The 
model reveals the respective roles in the different size range 
of two basic effects on nanoscale thermal transport, i.e., the 
intrinsic size effect is predominant at the size of about 5-100 
nm, and the surface scattering effect is dominant at the larger 
size. The theory agrees well with the existing experimental 
results for Si and GaAs nanowires and thin films. 
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